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The dc Josephson effect in a superconducting quantum point contact, where supercurrent flows 
through a small number of channels, is reviewed. The central role of Andreev levels is emphasized 
which carry the whole supercurrent in short symmetric Josephson junctions including tunnel junc- 
tions. A simple intuitive view of the dc Josephson effect in a quantum point contact is given in 
terms of multiple Andreev reflections. The quantization of the critical current in superconducting 
quantum point contacts is briefly discussed. 
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I. INTRODUCTION 

Recent surprising progress in micro-fabrication techniques has enabled us to make tiny Josephson junctions with 
artificial geometry. This opens a door to explore various mesoscopic effects in the transport properties of small 
superconducting junctions. In this paper we review the dc Josephson effect in a superconducting quantum point 
contact Jl|,||, a Josephson junction consisting of two superconductors coupled through a few quantum channels (Fig. 1). 
Such a superconducting quantum point contact can be created in a superconductor-2DEG-superconductor junction 
(2DEG: two-dimensional electron gas) j| as well as in a mechanically controllable break junction [0]. Main focus 
will be on current-carrying bound states that appear near a quantum point contact. Our discussion is based on the 
Bogoliubov-de Gennes equation, which is a two-component Schrodinger equation describing electron-like and hole-like 
quasiparticle excitations in superconducting systems. We show that the Bogoliubov-de Gennes equation has solutions 
describing bound states that are localized around the junction and decay exponentially into the superconductors. 
These bound states play a central role in the dc Josephson effect. The physical picture that emerges out of the 
solution is that the supercurrent flows as a result of multiple Andreev reflections |5| that transform an electron into a 
hole or vice versa. The bound states represent standing waves of electron and hole waves going back and forth between 
two superconductors. This view might look slightly different from a traditional picture of the dc Josephson effect 
explained in textbooks, according to which the supercurrent, calculated in lowest order in the tunneling Hamiltonian, 
flows as a result of tunneling of Cooper pairs. In fact, these views are simply two sides of a same coin. 
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FIG. 1. Superconducting quantum point contact. The supercurrent flows through a narrow constriction in the normal region. 

The importance of the bound states caused by the multiple Andreev reflections was first recognized in the early 
1970s in the studies of the dc Josephson effect in long superconductor-normal metal-superconductor junctions 
Until the 1990s, however, it was not fully appreciated that the presence of current-carrying bound states is a general 
feature in any symmetric Josephson junctions including weak links, superconducting quantum point contacts, and 
conventional tunnel junctions. The bound states in a short junction or a point contact, whose energy depends on the 
phase difference of the superconductors, were explicitly constructed by solving the Bogoliubov-de Gennes equation 
in 0,1). It is a curious fact that almost no attention was paid to the bound states in earlier studies of short 
superconducting weak links fll]-|l3|| , in which the bound states show up silently as poles of the Green's function on the 
real axis. In fact, the Andreev bound states play a central role for various situations in mesoscopic superconducting 
systems. 
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In this paper we will concentrate on the dc Josephson effect of s-wave superconductors. When a finite voltage is 
applied across the junction, the loop of the multiple Andreev scattering process is no longer closed, and the solutions 
of the Bogoliubov-de Gennes equation contain infinitely many terms each representing Andreev-reflected electron 
wave or hole wave 14-21 . For sufficiently small applied voltage, however, the current- voltage characteristics can 
be described in terms of the quasi-stationary bound states whose energy is changing slowly with time pcj . We will 
not discuss further the ac Josephson effect and the current-voltage characteristics of superconducting quantum point 
contacts in this paper. 

This paper is organized as follows. In Sec. II we introduce our basic equation, the Bogoliubov-de Gennes equation, 
and a continuity equation for electric charge. We solve a simplified model of a Josephson junction in Sec. Ill to 
illustrate the basic idea of the supercurrent carried by Andreev bound states. We point out that a midgap surface 
state in a ci-wave superconductor is essentially the same as the bound state in a short Josephson junction. We then 
explain a scattering approach to calculate the dc Josephson current in Sec. IV. 



II. BOGOLIUBOV-DE GENNES EQUATION 



We begin with the Bogoliubov-de Gennes equation [E2j 
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where U(r) represents both confining potential and scattering potential at interfaces of a superconductor and a normal 
metal, and p, is the chemical potential. This is the basic equation for mesoscopic superconducting systems. The upper 
component u(r) corresponds to electrons and the lower one v(r) represents holes. In principle the pair potential 
A(r) is determined from the selfconsistency condition or the gap equation, A(r) = A(4 f |(r)Wj p (r)), where A(> 0) is 
a coupling constant of electron-electron attraction, and ( ) is an average in the equilibrium state. The field operator 
^ a of electrons with spin a can be expanded with the solutions of the Bogoliubov-de Gennes equation: 
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and {a nr7 , a n / a >} = 0. Note 
That is why we sum over the positive 



where the creation and annihilation operators a nr7 and satisfy {aw, a> n > a '} — &n,n'&a,a 
that (— v* u* ) t is an eigenfunction with energy — E n if (u n v n )* satisfies Eq. 
energy states only and introduce the second term in Eq. (J2J) . 

Following Ref. [ p3[ , we introduce the equation of continuity of electric charge, which are satisfied by the solutions 
of the Bogoliubov-de Gennes equation. We define 
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It is clear that P e and J e are electric charge density and electric current, respectively. The continuity equation for the 
electric charge reads d t P e + divJ e + 5 = 0, where the source term S is given by S = (4e/^)Im[A(r)wj(r, £)\pj(r, t)]. 
When the gap equation is satisfied, (S) vanishes, and we recover the ordinary continuity equation of the charge for the 
averaged quantities (P e ) and (J e ), where ( ) is an average in equilibrium. When we calculate the Josephson current 
in a model where the pair potential satisfies the gap equation only approximately, we have to take into account the 
contribution from S The source term can be ignored when the current is calculated inside the normal region 

where A(r) = 0. 



III. CURRENT-CARRYING BOUND STATES 



We shall consider a simplest model of a Josephson junction that shows the essential features of the Josephson effect. 
It is a one-dimensional model where the left and right superconductors have the pair potentials of the same magnitude 
but with different phases: 
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We assume the normal region to be thin and consider an extreme case where its width is infinitesimal. In general we 
expect some scattering process to be present either inside the normal region or at the super-normal interfaces. We 
therefore introduce a scattering potential 



U{x) = VS(x) V>0 



(5) 



to describe such a scattering. The model with (Q) and (||) is the simplest non-trivial model for the Josephson effect. 
The strength of the (^-function potential is characterized by a dimensionless parameter Z = mV/h 2 kF, where fcj? is the 
Fermi wave vector. When 2< 1 the model describes a short superconductor-normal metal-superconductor junction. 
In the opposite limit Z>1, the model can be regarded as a toy model of a superconductor-insulator-superconductor 
tunnel junction. 

We are interested in a state in the energy gap, < E < Aq. In this energy range a wave function has to decay 
exponentially for \x\ — > oo. We can write the wave function as 
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where 



k h = k* e , u = Y^(l + ^), v = u* , T) = \J~A\ - E 2 . 



The coefficients a B , b B , c B , and d B are determined from the matching condition at x = 0: 

2mV 
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The eigenenergy is determined from the condition that Eq. (0) has a nontrivial solution. Within the approximation 
k e « kh ~ kp we find jl0[ 



E = E_ 

The solution of Eq. ([?]) with this energy is 
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Ncry'T+Z 2 v/cos 2 ^/ 2 ) + Z 2 ~ <JCOs(<p/2) , b B = iaNZ{l - iZ), 



c B = N(l-iZ) v/cos 2 ((p/2) + Z 2 -acos(p/2) , d B = iNZ^T+Z 2 , 
where <p — 9 B — 8l, N is the normalization constant 
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and a = sgn(y) for — ir < ip < n. The wave function is normalized such that [\u B (x)\ 2 + |t!s(a:)| 2 ] dx = 1. Once 

the wave function of the positive energy state is known, the eigenstate with energy — E B is obtained as i> B {x) = 
(— v B (x), u* B (x)) t . The field operator is then expanded as 
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Let us assume for a while that the states above or below the energy gap (\E\ > Aq) do not contribute to the 
supercurrent. We will see later that this is indeed the case for the model we are dealing with. We can calculate the 
supercurrent as the average of J e at x — +0 in the equilibrium state at temperature T: 
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where /(-E) = [exp(E/k B T) + 1] 1 . Since the transmission probability through the J-function barrier Eq. (|5|) is 
T = 1/(1 + Z 2 ), the normal-state conductance of this model is Gn = 2e 2 /[h(l + Z 2 )]. We can thus rewrite Eqs. (||) 
and @ as 
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When T < 1 Eq. (lit) reduces to 



I 



-Gn n ° sin(y;) tanh [ 



2c 



2fc B T 



(12) 



which is the well-known Ambegaokar-Baratoff formula |25[|. In the opposite limit T = 1, Eq. (lib) becomes 

/ = - (TV v-^ sin 1^/2 jt anh[ ms( ,:/2 ) | . (13) 



e ^' ' \2k B T 

which coincides with the result obtained by Kulik and Omel'yanchuk |]IT| for a point contact Josephson junction. In 
this limit the zero-temperature Josephson current becomes maximum at tp = it, not at ip = tt/2 like in the tunnel 
junctions. The bound state energy is simply given by E B — Ao cos(</?/2), which equals the Fermi energy E B — at 

(f = 7T. 



Equation (lit) thus interpolates between Eq. (|12|) of tunnel junctions with very small transmission probability 
and Eq. ( |l3| ) of point contacts with perfect transmission. The same expression of the supercurrent interpolating 



the two limits has been found by many people [g(|,|l2|,[f3 27 1. The phase-dependent bound state energy (|3|) and its 
implication to the Josephson effect in short junctions were emphasized in ]To|| . The results were then generalized to 
the multichannnel contacts in pgj ]. 

Equation (JsJ) shows that there is always a bound state (Andreev level) with E B > unless the phase difference 
tp is zero. This is true even for the tunnel junctions in which T <^ 1 and Ao — E B <C Ao. In the above model two 
bound states (E = ±E B ) appear in the energy gap \E\ < Aq but disappear when tp = 0. This is a general property 
of short s-wave junctions where hvp/ L ^> Ao with L being the length of the normal region (L = in our model). In 
the opposite regime where Tivf/L Ao, many (~ hvp/ AqL) Andreev levels exist in the gap for any tp Iq-pJ]. These 
Andreev levels are easily understood as discrete levels of normal electrons confined in a quantum well of offdiagonal 
pair potential. They are similar to the discrete energy levels in a vortex in an s-wave type II superconductor | p9[ . A 
nontrivial observation we made in the above calculation is that, when ip ^ 0, two Andreev levels survive even in the 
limit L — > 0, in sharp contrast with a case where a particle is confined in a quantum well of diagonal potential barrier, 
in which no bound state exists in the limit L —> 0. 
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FIG. 2. Current-carrying bound state. The bound state is a manifestation of the multiple Andreev reflections and carries the 
whole supercurrent in a short Josephson junction. The filled (open) circle represent an electron (a hole). The ovals containing 
two circles are Cooper pairs. 
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The bound state iPb{x) is formed as a result of resonant multiple Andreev reflection process. Equation ([lO|) shows 
that this state carries nonvanishing supercurrent. This can be physically understood as follows (Fig. ||). Each time 
an electron is Andreev reflected into a hole, a Cooper pair is effectively generated. Therefore the bound state, which 
represents an infinite loop of Andreev reflections (electron— > hole— >electron- • ■), serves as a pump that transfers Cooper 
pairs from one superconductor to the other. When T rs 1 the bound state energy has a strong dependence on the 
phase difference ip, which can be controlled externally by the current source. It would be very interesting if these 
current-carrying Andreev levels are observed experimentally using, for example, the scanning tunneling spectroscopy. 

Finally we comment on midgap states that appear near the surface of a d-wave superconductor when the surface 
is not parallel to a crystal axis (Fig. ||) [Q . In fact the mechanism of the formation of the zero-energy midgap state 
is the same as the Andreev bound states we discussed above. As shown in Fig. pi an electron traveling towards a 
surface is reflected back into the d-wave superconductor and is subsequently Andreev-reflected into a hole by the 
positive pair potential. In the next step the hole follows the same path backwards, reflected at the surface, and finally 
Andreev-reflected into another electron by the negative pair potential. The analogy to the Josephson junction is now 
obvious. The surface of the d-wave superconductor plays a role of the poi nt co ntact with T = 1, and the sign change 



in the pair potential corresponds to the phase difference tp = tt. Equation (11a) tells that there must be bound states 
with Eb = localized near the surface, which are the midgap states found in Rcf. |50|. They give rise to a peak in 
differential conductance of a normal-insulator-d-wave tunnel junction |3l[| . The Andreev levels also play an important 
role in a Josephson junction made of an s-wave superconductor and a d-wave superconductor |5^j53|. It is known 
that in such a junction there is a chance that Andreev levels appear even in the lowest-energy state, thereby leading 
to a ground state with broken time-reversal symmetry |p4f . 
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FIG. 3. Midgap state formed near the surface of a d-wave superconductor. 



IV. JOSEPHSON CURRENT FROM SCATTERING AMPLITUDES 



We shall introduce a general expression of the dc Josephson current under the assumption that we can ignore the 
spatial variation of the pair potential. Here we assume 

A(r) = /^ eiS ; * <0 < (14) 
V ' \& R e l9R x>L. V ' 

In the normal region (0 < x < L) the order parameter vanishes. We also assume that the potential U(r) is independent 
of x in the superconductors. We emphasize that this step-function pair potential ( |T4| ) is a good approximation for 
mesoscopic Josephson junctions of our concern, where bulk superconductors are coupled via a few conduction channels. 
The use of the step-function model is justified by the following reasons. First, the presence of a point contact does 
not seriously affect the order parameter of bulk superconductors, because the width of the constrictions is much 
smaller than the superconducting coherence length £, which is the characteristic length scale for the variation of the 
order parameter. Second, if the normal region is made of a material different from superconducting electrodes, the 
contact resistance of the SN interfaces is not small due to the Fermi- velocity difference and/or due to the presence 
of a Schottky barrier. In this case also the pair potential is not reduced significantly near the point contact. More 
generally, when the resistance of the junction is much larger than those of superconductors, we may safely use the 
step- function model |55|. 
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Since the dc Josephson current is an equilibrium current, it is given by the thermal average of J e . We calculate it 
at x = to avoid complication due to the source term S: 



(15) 



2=0 



where G w (r,r') is the temperature Green's function at frequency u) n = 7rfcsT(2n + 1) in the Nambu space (2x2 
matrix) [BfJ. Our assumption that the order parameter A(r) and the potential U(r) have no dependence on x in the 
superconductors, is strong enough that we can immediately write down G u , n (r ) r') at x < in terms of scattering 
amplitudes. Its explicit expression can be found in Ref. p4| . After substituting it into Eq. Jl5|), we obtain 



— ^[ajjfaiun) - ajj(-<p,iuj n )], (16) 



h(3 ^ n ni 

where Q n L — \J w% + A L and we have used the Andreev approximation where the velocities of quasiparticles are set 
equal to the Fermi velocity. 

The quantity a,jk(ip,iw n ) is a scattering amplitude for the process in which an electron-like quasiparticle of jth 
channel traveling from the left to the junction is reflected back as a hole-like quasiparticle of fcth channel . Equation 
( |l6| ) expresses the supercurrent in terms of the scattering data cijj, just like the Landauer formula does the normal- 
state conductance. An important difference is that the supercurrent is proportional to the scattering amplitude, not to 
the probability. Equation (|l^) is clearly an odd function of tp, and can be expanded formally as / = Ylm= l ^ m si n (m<^). 
The term proportional to sin(mt/?) corresponds to mth order term in the tunneling Hamiltonian approach that is a 
contribution from the simultaneous tunneling of m Cooper pairs. Equation ( |l6| ) is thus a useful, nonperturbative 
formula for the dc Josephson current |$7j . 

Let us apply it to the simple one-dimensional model we analyzed in the last section. Solving the Bogoliubov-de 
Gennes equation for E > Aq with potentials given by (||) and (||) gives 

^ _ A [E sin 2 (y/2) + ito cos(y>/2)] 
' E 2 (l + Z 2 )-A 2 [cos 2 (<p/2) + Z 2 ] • 1 ' 



Note that the scattering amplitude has poles at E = ±Er- Substituting (jlTj) with E = iuj n into (|16[), we then find 

_eA| v sinjp 

HP j£ w 2 (l + Z 2 ) + A 2 [cos 2 (^/2) + Z 2 ]' 1 ' 

which agrees with Eq. (|l0|) after the frequency sum is evaluated by reducing into a contour integral. In this calculation 
all the contributions to the integral come from the residues at E = ±Er- This proves that the Andreev bound states 
carry the whole supercurrent in short symmetric Josephson junctions. We note that this is not true for long weak 
links where hvp/L <c Ao and for asymmetric junctions where A^ < A^. Let us examine the latter case in 
more detail. We consider a one-dimensional model where the pair potential is given by Eq. ( |l4| ) with L = and the 
diagonal potential is (||). A solution of the Bogoliubov-de Gennes equation for this model gives a scattering amplitude 
for E > A R > A L , 

E(Ar cos ip — Al) — iAr^Il sin ip 
^ E) = £ 2 + W^-A fl A L cos„ ' ^ 

where VLr {l) = (E 2 - A R{L) )^ 2 and K = 1 + 2Z 2 . From Eq. (||) we obtain 

I=-—\^ ArAl sin <p 

hp wi + Kn nR a nL + ArA l cos 

where 0„jj(£) = (uj 2 + A 2 ^^) 1 / 2 . This result agrees with the Ambegaokar-Baratoff formula for a tunnel junction 
in the limit K ~3> 1 and with the Kulik-Omel'yanchuk formula [ |TT| | for a point contact with perfect transmission in 
the limit K — > 1. The summation can be reduced to an integral over A^ < \E\ < Ar plus contributions from residues 
at energies of bound states at \E\ < Al. The integral over A < \E\ < A^ corresponds to the process shown in Fig. || 
and carry most of the supercurrent. 
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FIG. 4. Scattering states carrying supercurrent in a short asymmetric junction. 



V. QUANTIZATION OF SUPERCURRENT 

In a superconducting quantum point contact, the critical current or the maximum value of the dc Josephson current 
can take a quantized value in units of eAo/h. This was first predicted theoretically and confirmed by experiments 
HQ. The effect can be understood easily by the following argument. The critical current of a Josephson junction is 
generally proportional to its normal-state conductance, as seen in Eq. ( |llb| ). It is well known that the conductance of 
a quantum point contact fabricated in a 2DEG is quantized in units of 2e 2 /h |^,^|: The conductance is simply 2e 2 /h 
times the number of open channels at the constriction. We can thus naively expect that the critical current should 
be quantized as well. For a more quantitative discussion we can use the adiabatic approximation where scattering 
between channels is neglected. The transmission probability through the quantum point contact is then 7j = 1 for 
open channels and T 3 ■ = for closed chan nels. For a short symmetric superconducting quantum point contact where 
Hvf/L ^> A , we obtain, from Eqs. (lib) and (|l^), 



/ = S m(<p/2) tanh cosfe>/2)) , (21) 

where N is the number of open channels at the point contact. We thus find that the critical current is a multiple of 
eAo/h at zero temperature. The number of open channels is determined by the width of the constriction, which can 
be controlled externally, for example, by changing a gate voltage for a 2DEG. This was done in an experiment by 
Takayanagi et al. jj). 

Equation (|2l| ) is valid only for short junctions. In the opposite case where Tivp / L <C Ao, the unit of the quantization 
is different. The supercurrent for a long superconducting quantum point contact is given by 



N 

e 



( 22 ) 
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for \<p\ < 7T at zero temperature, where Vj is the longitudinal component of the Fermi velocity in jth channel. In 
this case the critical current is not exactly quantized in contrast with the short-junction limit. This is nothing but a 
quantum limit of long SNS junctions obtained in early 1970s [Mil. 



VI. CONCLUSIONS 



In this paper we have reviewed the dc Josephson effect in mesoscopic Josephson junctions. Our discussion is based 
on a simple analysis of the Bogoliubov-de Gennes equation. In short junctions including superconducting quantum 
point contacts the whole supercurrent is carried by Andreev bound states in the energy gap that appear as a result 
of multiple Andreev reflections. The wave function of a bound state is constructed explicitly for a simple (S-function 
model by solving the Bogoliubov-de Gennes equation. The bound state energies are sensitive to the difference of the 
phases of superconductors, and this sensitivity is the origin of the Josephson effect. Our simple physical view of the 
Josephson effect in terms of the current-carrying bound states is very general and applicable to many situations. 
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